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TIRE PRESSURE ESTIMATION 



FIELD OP THE INVENTION 

The present invention relates in general to a method of, a sys- 
tem for and a computer program product for estimating the oc- 
currence of a specific tire pressure deviation between actual 
and nominal pressure values for one or a plurality of wheels. 

BACKGROUND OF THE INVENTION 

Driving safety and- driving economy play an important role in 
automobile design. Tire pressure monitoring systems, which con- 
tribute to safe driving and enhanced fuel efficiency, are in- 
creasingly integrated in automobile equipment. Two different 
techniques, namely direct and indirect tire pressure measure- 
ments, are distinguished in the prior art. Direct measurements, 
which are very precise, require comparatively expensive equip- 
ment integrated within the tire. Indirect measurement systems, 
which may be based on anti- locking systems (ABS) , are easily 
implemented and cost-efficient but up to now less precise than 
direct measurement techniques. 

Indirect measurement techniques comprise wheel radius analysis 
(WRA) and wheel vibration analysis (WVA) . 

In wheel radius analysis, the wheel radii are estimated by com- 
paring the wheel speeds of each vehicle wheel (which may be 
measured by wheel speed sensors of the ABS) with the vehicle's 
actual driving velocity. Constant driving velocity provided, 
the wheel speed increases when the wheel radius decreases due 

* 

to tire pressure loss . 

In wheel vibration analysis, wheel vibration phenomena (which 
may be measured by wheel speed sensors of the ABS) are detected 
within the time dependent behavior of the rotational velocity 
of an individual vehicle wheel. Tires in contact with the road 
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surface during driving are exposed to vibrations. Here, the 
fact is exploited that the spectral properties (e.g. the energy 
distribution in different frequency bands) of a wheel angular 
velocity signal depends on the tire pressure of the individual 
vehicle wheel . 

Traditional tire pressure estimating systems rely on either 
wheel radius analysis or wheel vibration analysis. Both systems 
have particular advantages and deficiencies. Estimating tire 
pressures with wheel radius analysis is easily implemented but 
affected by vehicle load changes or driving forces during cor- 
nering, acceleration and braking, which also lead to wheel ra- 
dius changes. Wheel vibration analysis, which requires in- 
creased computing performance, is very precise in a low and me- 
dium velocity range, but shows deficiencies in a high velocity 
range of the driving velocity. 

JP 5-133831 discloses a WVA which identifies a pronounced low 
frequency resonance which can be measured when the vehicle ve- 
locity is low. Typically, the. corresponding resonance frequency 
lies in the range of about 30 to 50 Hz. When the vehicle veloc- 
ity increases to a high velocity range, the low resonance peak 
looses intensity. In this velocity range the JP 7-21723 and EP 
0925960 further identify a high frequency resonance which is 
suitable for the tire pressure estimation. The frequency of 
this resonance is in the range of about 60 to 90 Hz. However, 
for too high vehicle velocities, it becomes increasingly diffi- 
cult to use vibrational analysis for tire pressure estimation. 
FIGS. 12A and B show two power spectra of the rotational ve- 
locities which correspond to a low and a high vehicle speed, 
respectively. The low (FIG. 12a) and high (FIG. 12b) resonance 
peaks can easily be identified. 

The low resonance peak is generally ascribed to a torsional 
resonance in the rotational direction of a tire. When the tire 
air pressure drops the spring constant in the torsional direc- 
tion decreases. The high resonance peak is either ascribed to 
secondary components of the torsional resonance frequency or to 
tread pattern effects. FIG. 13 shows a typical relation of the 
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tire pressure as a function of the resonance frequency (high 
resonance peak) . 

EP 0 925 960 and WO 01/87647 disclose tire air pressure esti- 
mating systems implementing both analysis techniques, WVA and 
WRA, to compensate for the individual deficiencies of these 
different analysis techniques. Hereby, the system described in 
EP 0 925 960 switches between the two analysis techniques de- 
pending on the current driving parameters and conditions 
(vehicle velocity, resonance signal intensity, vibration input 
intensity, brake switch signal, etc.). When the vehicle is in a 
low or medium velocity range, where WVA is efficient, the tire 
air pressure estimating is based on the output of the WVA. In 
the high velocity range, where WVA is less efficient, it relies 
on the output of the WRA. The tire pressure estimating system 
disclosed in WO 01/87647 derives for each vehicle wheel two in- 
dividual confidence level values corresponding to the outputs 
of the two analysis techniques, WVA and WRA, respectively. If 
the outputs of the WVA and the WRA correspond to each other, 
then the sum of the corresponding two individual confidence 
values is compared with a first threshold value. Otherwise, the 
higher value among the two individual confidence level values 
is chosen and compared with a corresponding threshold value 
(one threshold value for the WVA and WRA, respectively) . A tire 
pressure indication signal indicating a tire pressure loss is 
generated when one of the three threshold conditions is ful- 
filled. 

OBJECT OF THE INVENTION 

The object of the invention is to improve the liability of the 
estimation of the occurrence of a specific tire pressure devia- 
tion for a wheel . 

SUMMARY OF THE INVENTION 

According to a first aspect the invention achieves this object 
by a method of estimating the occurrence of a specific tire 
pressure deviation between actual and nominal pressure values 
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for one or a plurality of wheels. The method comprises subse- 
quently obtaining one or more wheel radius analysis measurement 
values from a wheel radius analysis component, wherein the 
wheel radius analysis measurement values are related to single 
wheel radius values of which each is indicative of the wheel 
radius of a particular wheel. It further comprises subsequently 
obtaining one or more wheel vibration data values from a wheel 
vibration analysis component, wherein each of the wheel vibra- 
tion data values is indicative of a vibration phenomena in the 
time dependent behavior of the rotational velocity of a par- 
ticular wheel. Finally, it comprises calculating one or more 
tire pressure output values on the basis of both the wheel ra- 
dius analysis measurement values and the wheel vibration data 
values wherein each tire pressure output value is indicative of 
the specific tire pressure deviation for a particular wheel. 

According to a second aspect the invention achieves this object 
by a system for estimating the occurrence of a specific tire 
pressure deviation between actual and nominal pressure values 
for one or a plurality of wheels. The system comprises a compo- 
nent for subsequently obtaining one or more wheel radius analy- 
sis measurement values from a wheel radius analysis component, 
wherein the wheel radius analysis measurement values are re- 
lated to single wheel radius values of which each is indicative 
of the wheel radius of a particular wheel. It further comprises 
a component for subsequently obtaining one or more wheel vibra- 
tion data values from a wheel vibration analysis component, 
wherein each of the wheel vibration data values is indicative 
of a vibration phenomena in the time dependent behavior of the 
rotational velocity of a particular wheel. Finally, it com- 
prises a component for calculating one or more tire pressure 
output values on the basis of both the wheel radius analysis 
measurement values and the wheel vibration data values wherein 
each tire pressure output value- is indicative of the specific 
tire pressure deviation for a particular wheel. 

According to a third aspect the invention achieves this object 
by a computer program product including program code for carry- 
ing out a digital signal processing method, when executed on a 



WO 2005/005173 



PCT/EP2003/007283 



- 5 - 



computer system, for estimating the occurrence of a specific 
tire pressure deviation between actual and nominal pressure 
values of one or a plurality of wheels. The method comprises 
subsequently obtaining one or more wheel radius analysis meas- 
urement values from a wheel radius analysis component, wherein 
the wheel radius analysis measurement values are related to 
single wheel radius values of which each is indicative of the 
wheel radius of a particular wheel. It further comprises subse- 
quently obtaining one or more wheel vibration data values from 
a wheel vibration analysis component, wherein each of the wheel 
vibration data values is indicative of a vibration phenomena in 
the time dependent behavior of the rotational velocity of a 
particular wheel. Finally, it comprises calculating one or more 
tire pressure output values on the basis of both the wheel ra- 
dius analysis measurement values and the wheel vibration data 
values wherein each tire pressure output value is indicative of 
the specific tire pressure deviation for a particular wheel. 

Other features inherent in the tire pressure estimating method 
and system and the corresponding computer program product are 
disclosed or will become apparent to those skilled in the art 
from the following detailed description of embodiments and its 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the accompanying drawings, 

FIG. 1 shows a schematic diagram of a first embodiment 



of the invention; 



FIG. 2A, B 



show two contour plots of the pressure indica- 
tion value ^ for two different sets of tuning 
parameters; 



FIG. 3 



shows a flow diagram of the process steps of an 
implementation of the first embodiment; 



FIG. 4 



shows a flow diagram of the process steps of a 
preferred tire radius transformation which is 
performed in the first embodiment; 
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FIG. 5 shows a schematic diagram of a second embodi- 

ment of the invention; 

FIG. 6 shows a flow diagram of the process steps of an 

implementation of the second embodiment; 

FIG. 7 shows a schematic diagram of a third embodiment 

o f the invent i on ; 

FIG. 8 shows a flow diagram of the process steps of an 

implementation of the third embodiment; 

FIG. 9 shows a schematic diagram of a fourth embodi- 

ment of the invention; 

FIG. 10 shows a schematic diagram of a load balance 

compensation used in the fourth embodiment; 

FIG. 11 shows a flow diagram of the process steps of an 

implementation of the fourth embodiment; 

FIGS. 12A, B show two power spectra of the rotational ve- 
locities which correspond to a low and a high 
vehicle speed, respectively, according to the 
prior art; 

FIG. 13 shows a typical relation of the tire pressure 

as a function of the resonance frequency ac- 
cording to the prior art; and 

FIG. 14 shows four different graphs, one for each of 

the four tires, in which ( A/) , Ar f ) parameter 
pairs obtained for each single tire i are plot- 
ted during a test drive under normal single 
tire conditions. 



DETAILED DESCRIPTION OF THE EMBODIMENTS 



Below, rotational speed sensors as used for measuring rota- 
tional velocities of vehicle wheels, a wheel radius analysis 
(WRA) component and a wheel vibration analysis (WVA) component 
are explained in more detail which are common parts of the 
various embodiments of the invention. 
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[Rotational speed sensor] 

Rotational speed sensors are often used for measuring the angu- 
lar motion of a rotating axle. The wheel axles in vehicles, the 
cam shaft in engines and the motor axles in robots are some ex- 
amples for the application of these sensors. Such a sensor may 
comprise a transducer (for example an inductive or optical 
transducer, a Hall sensor, etc.) cooperating with the teeth of 
a toothed wheel. 

Every time a tooth passes the transducer the latter one gener- 
ates a trigger signal. The rotational speed is determined by 
measuring the time elapsed between two adjacent trigger sig- 
nals. Between two adjacent trigger signals the toothed wheel 
rotates about the angle a^lxIN, where N is the number of 
teeth of the wheel. The trigger signals may be sent to a micro- 
processor unit where they are converted to angular velocity 
values. Often rotational speed sensors built in the anti- 
locking system (ABS) of a vehicle are used to provide these an- 
gular velocity values. 

A more detailed description of rotational speed sensors and 
their application in the context of tire pressure estimating 
systems can be found in EP 0925960. 

For each of the vehicles wheels, there is provided a permanent 
input of corresponding angular speed values to a wheel radius 
analysis component and to the wheel vibration analysis compo- 
nent. The input data values may be provided as time domain sam- 
pled data sets ( equidistant ly spaced in time) or as event do- 
main sampled data sets (non-equidistantly spaced in time) . 

[Wheel radius analysis (WRA) component] 

A WRA is based on the fact that the wheel speed of a wheel de- 
pends on the respective wheel radius. The smaller the wheel ra- 
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dius is, the larger is the corresponding wheel speed. A WRA 
component may in particular provide differences of wheel speeds 
from pairs of vehicle wheels. The wheel speeds of each wheel 
are thereby measured by the rotational speed sensors of each 
wheel. In general, the WRA component outputs a number of WRA 
measurement values AR of which each is related to either one, 
two, three or all four single tire radius values Ar i indicating 
the changes in the wheel radius of a particular wheel i. 

According to one preferred embodiment, the WRA component pro- 
vides a four-dimensional vector 

AR = [A„,A XL ,A 1H ,A ni ] r . da) 

where A^A^A^ and are the differences in the wheel radii 

of the front left and front right, front left and rear left, 
front right and rear right, rear left and rear right wheels, 
respectively. The vector AR is related to the single tire ra- 
dius values Ar i via the following equation: 



-—HAr , (2a) 

where r 0 is a known nominal wheel radius and each of the indi- 
ces 1,2,3,4 of Ar f corresponds to a respective vehicle wheel: 1 
= front left (FL) , 2 = front right (FR) , 3 = rear left (RL) , 4 
= rear right (RR) . 

It is noted, that the parameters , A^ , A YF and A^ in fact 
depend on the tire pressures of the corresponding wheels, but 
A„ and A w may also reflect vehicle load changes or surface 

XL -*R 

changes and similarly, and A ra may react on driving forces 

(acceleration, braking, forces in curves, etc.). 



-11 0 0 
-10 10 
0-101 
0 0-11 



Ar, 

>4. 



WO 2005/005173 



PCT/EP2003/007283 



9 - 



If the vehicle does not include any gyro which outputs yaw 
rates of the driving vehicle an alternative WRA component may 
provide a three-dimensional vector 



Afl^A^A^-A^A^f 



(lb) 



In this case, the matrix H of equation (2a) may be adopted as 

follows : 



-11 1-1 
-1 0 1 0 
0-101 



(2b) 



In general terms, the WRA component may consider any arbitrary 
function of the wheel radius values. For example, the WRA com- 
ponent applies the relationship 

R = f(r u r 29 r S9 r 4 ) 
(lc) 

where f is an one, two, three, or four-dimensional function 
that relates the wheel radius values to an observable (vector- 
valued) quantity jR. Applying a first order Taylor series ap- 
proximation, the expression R = f(r u r 29 r 3 ,r 4 ) can for small wheel 
radius deviations be written as 



R*f(r 0> r 0 ,r 07 r Q ) + 



(Id) 
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The gradient matrix 

M_ sL <L Ml 

dr x dr 2 9r 3 9r 4 
(2d) 

should be evaluated using the nominal value r 0 for the wheel 
radius. Consider for example the following quantity: 

R = H-?±. (le) 
Application of equation (Id) then shows that 




(If) 



The outputs of a WRA component are not limited to the exempli- 
fied outputs of equations (la-f ) . For each specific WRA compo- 
nent, the Taylor series approach (equation (Id)) can be applied 
to construct a corresponding gradient matrix H which relates 
the WRA measurement values AR to the single tire radius values 
A/;. . 

[Wheel vibration analysis (WVA) component] 

The WVA component used in the embodiments of the present inven- 
tion can be of any type contemplated, in particular of those 
types as mentioned in the introductory part of the application. 
For a brief summary of the principles of WVA reference is made 
as well to the prior art as discussed above. 

The WVA component outputs wheel vibration data values (for ex- 
ample the energy distribution in different frequency bands in 
the time dependent behavior of the rotational velocity of a 
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particular wheel (/ ) , one or more resonance frequencies of a 
specific wheel (z ) or any other suitable measure) . In the fol- 
lowing, these WVA component outputs are generally denoted as 
vibration phenomena values. 

The WVA component used in the embodiments of the present inven- 
tion provides a four -dimensional vector Af whose components 

indicate the deviation between a currently measured vibration 
phenomena value and a corresponding calibration value: 



The indices i = 1, 2, 3 and 4 here again denote the respective 
wheels FL, FR, RL and RR (see above) . A tire pressure drop in a 
particular tire will lead to a decreasing value of the cur- 
rently measured vibration phenomena value. The corresponding 
component of Af will hence become significantly negative when 
the tire pressure decreases. 

[First embodiment] 

FIG. 1 shows a schematic diagram of a first embodiment accord- 
ing to the invention. In a preprocessing step 106 (according to 
equation (2a), see above), the values AR i ( j = YF, XL, XR, YR ) ob- 
tained from a WRA 104 are transformed to single tire radius 
values Ar { (i = FL,FR 9 RL,RR) . Then, probabilities Pf and f> r are 
computed from the vibration phenomena values isf. obtained from 
a WVA 102 in step 108 and from the single tire radius values 
Ar t in step 112, respectively, on the basis of a Gaussian cumu- 
lative probability distribution function by treating the values 
A/] and Ar { as independent Gaussian distributed random vari- 
ables. In a further step 110, a weight factor is calculated 
to penalize disagreement between the outputs Af f and AR. from 
the WVA and WRA. Based on the obtained values Pf , P { r and W/ r , 
a tire pressure indication value 77, is computed in step 114 , 
which is then compared with a predetermined threshold value in 




(3) 
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step 116. The above steps 106 to 116 are described in more de- 
tail in the following. 

In step 106, the measurement vector AR is preferably trans- 
formed into the vector Ar , whose elements describe the changes 
in wheel radius of each particular vehicle wheel (see also 
above) . This transformation involves the inversion of the ma- 
trix given in equation (2a) , and this inversion is non-trivial 
since said matrix is singular. The matrix is singular since the 
WRA component is not able to detect pressure deviations that 
simultaneously occur in all four wheels. Assuming that the 
wheel radius cannot increase, a pseudo- inversion of equation 
(2a) is however obtained by the following transformation: 

r^vtft.AK , (5) 
Ar = r tmp -mzx(r tmp )\\ 1 1 if . (6) 



Here, the vector r is only used in an intermediate step. The 

* imp 

matrix is the pseudo- inverse of the matrix H and is de- 

fined as follows: 



-3 -3 -1 -1 

3 -1-3 1 

-13 1-3 

113 3 



(7) 



In steps 10 8 to 114 the following computations are performed, 
in order to estimate the probability (or confidence level) 7jj 
(a possible pressure output value of the tire pressure estima- 
tion system) for the hypothesis "no deflation of tire i". 

In step .108, the probability Pf of a vibration phenomena de- 
viation for tire i which is larger than Af { is computed from 



(8) 
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Similarly in step 112, the probability P f r of a .wheel radius 
deviation for tire i which is larger than At; is computed from 



(9) 



Here, the function <D(x) is the cumulative probability distribu- 
tion function (CDF) for a Gaussian random variable with stan- 
dard deviation a = 1 : 



(10) 



In equations (8) and (9) , a f and a r are the corresponding 

standard deviations of the vibration phenomena value and tire 
radius measurements, respectively- Preferably, the standard de- 
viations may be predetermined tuning parameters obtained from 
test drives . These tuning parameters are stored in a tuning pa- 
rameter memory 10 0. They stay constant during operation of the 
tire pressure estimating system. For those skilled in the art, 
various alternatives for choosing or measuring the standard de- 
viation parameters are easily contemplated, including those 
with dynamically adapted values. 

In step 110, an additional weight factor Wf is calculated in 
order to penalize disagreement between the outputs A/) and AR. 
from the WVA 102 and WRA 104, which is defined as follows: 



Wf r = exp 











f 


A/ ( Ar,ri 








•exp 




cr f a r \) 



(11) 



and whose effect will be explained in more detail below, cr, and 
cr 2 are further tuning parameters (cp. above) . 
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In step 114, the pressure indication value tj. (namely, the test 
quantity in the above hypothesis test) for a flat tire is com- 
puted as follows: 

5 tj. ^C-Pf-Pf-W? . (12) 

The constant C is another tuning parameter which is predefined 
and stored in the tuning parameter memory 100. It is chosen as 
.such that if there are no deviations from the calibration val- 
10 ues ( Af. = 0 , Ar. = 0 ) , the confidence level for the hypothesis "no 
deflation of tire i" is tj. = 1 . Negative values of the deviation 
parameters Af. and Ar { lead to a decrease of the pressure indi- 
cation value ' tj . The value ij i is an indicator for the statisti- 
cal significance of the occurrence of a specific pressure de- 
ls viation (e.g. flat tire condition, pressure drop about a given 
amount, pressure increase about a given amount) . 

In step 116, the pressure indication value tj. is compared with 
a critical threshold value rj c . If ij falls below rj c (*Jl<*7 c )# a 

20 tire pressure loss signal (alarm) is triggered. The tire pres- 
sure loss signal may for example be issued if rj i < 0.1. To 
avoid toggling between alarm/no alarm, the alarm is not reset 
if Tj f increases just slightly above the threshold value tj c - In- 
stead it is required that 77, increases to a value well above 

25 the threshold value rj c . The desired specific pressure deviation 
whose occurrence shall be monitored may be adjusted by the 
standard deviation parameters a f and a r as well as the criti- 
cal threshold value tj 0 . 

30 The effect of introducing the weight factor Wf is to enhance 
the decrease of 77, for substantially conforming deviation val- 
ues from the WVA 102 and WRA 104 ( Af/a f « Ar, ja r ) and to suppress 
the decrease of 77, for substantially non- conforming deviation 
values ( Af t /a f * Arja r ) . The strength of the effect can be tuned 
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with the parameters <j x and cr 2 / which preferably have prede- 
fined constant values obtained from test drive situations, 
similar to the parameters a f and a r . a x and <x 2 are adjusted as 

such that the desired triggering of the tire pressure loss sig- 
nal is obtained. 

To illustrate this effect, contour plots of the pressure indi- 
cation value ?]. for two different sets of tuning parameters a x 
and <j 2 are shown in FIGS. 2A and B. For the plot shown in FIG. 
2A, the parameters were chosen to be small ( cr, = 0 , cr 2 = 0.4 ) in 
order to introduce merely a small penalizing effect. In this 
case, the gradients in direction of the axes (distance of the 
contour lines near the axes) , where WVA and WRA disagree with 
each other, is similar to the gradient along the diagonal 
btfijijf- Arja r , where WVA and WRA agree with each other. For the 
plot shown in FIG. 2B, the parameters were chosen to be large 
(o-j = 0.5 / <r 2 = 1.0) in order to introduce a strong penalizing ef- 
fect. This' renders the gradient along the diagonal steeper than 
the gradient in direction of the axes. 

FIG. 3 shows a flow diagram of the method for calculating the 
pressure indication value ij . In a first step 3 00, the tuning 
parameters <j f , a r , a } , cr 2 and C are received from the tuning 
parameter memory 100. In a second step 302, the current vibra- 
tion phenomena deviation value A/) is received from a WVA com- 
ponent, and, in a third step 304, the current tire radius de- 
viation value ARt is received from a WRA component. In a fur- 
ther preprocessing step 306, the tire radius deviation value 
Afl, is transformed to a single tire radius deviation value Ar. . 
In two subsequent steps 308 and 310, the Gaussian CDF values 
Pf and P f r are computed from the quantities Af. , cr f and Ar. , a r , 

respectively. In a next step 312, the weight factor W/ r is com- 
puted from the quantities Af t , Ar i , a f , a r , a x and cr 2 - Finally, 
the pressure indication value ^ is computed from C, Pf , P[ 
and W/ r in step 314, which is compared with the predetermined 



WO 2005/005173 PCT/EP2003/007283 

- 16- 

threshold value tj c in step 316. If the pressure indication 
value 7j. is less than the threshold value rj c , a tire pressure 
loss signal is triggered in step 318. The procedure is repeated 
for every new data received from the WVA and WRA components. 

5 

FIG. 4 shows a flow diagram for the preferred preprocessing of 
wheel radius deviation values A/^ (step 3 06 in FIG. 3) . In a 
first step 400, the wheel radius deviation values A/?, are re- 
ceived from the WRA component. In a second and third step 4 02 
10 and 404, the predetermined pseudo- inverse matrix and the 

known nominal wheel radius r 0 / respectively, are read from a 
persistent memory (e.g. the parameter tuning memory 100 of FIG. 
1) . In step 406, a temporary radius deviation vector r^ p is 

calculated by computing the product Tq^H^-AR. In a subsequent 
15 step 408, the maximum element of the vector r is determined 

imp 

and then, in step 410, the single tire radius deviation values 
Ar f are computed therefrom (which are further processed in step 
310 of FIG. 3) . 



20 [Second embodiment] 



25 



FIG. 5 shows a schematic diagram of a second embodiment accord- 
ing to the invention. Here, the model underlying the tire esti- 
mation is based on the assumption that the wheel radius value 
A/^ and the wheel vibration value A/) depend on the same pres- 
sure parameters Ap, instead of treating these measurement val- 
ues as independent quantities. The tire pressure vector 



30 



Po.fl Pfl Pojtr Pfr Po,rl Prl Po,bjr Prr 

Po t FL 



IT 



PojTR 



PojlL 



PojiR 



(13) 



consists of four pressure deviation values Ap. , one for each 
wheel (i=FL,FR,RL,RR) . p i is the current pressure value of the 
wheel i to be determined, and p 0{ is the corresponding calibra- 
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tion pressure value. (The Ap t are possible pressure output val- 
ues of the tire pressure estimation system.) 

In a load change transformation step 504, the wheel radius val- 
ues AR i obtained from the WRA 502 are transformed into modified 

wheel radius values AQ in order to reduce the impact of vehi- 
cle load changes. Hereby, the four components of the vector 
AR are reduced to three components of the vector AR , which 

are defined as AR = [A^ A^-A^ A^] 7 (cp. equation (lc) ) . The 
quantity A^-A^ of the vector AR is considerably less sensi- 
ble to load changes than the quantities A^ and A^ of the 
vector AR . 

According to the model of this embodiment, the relation between 
on the one hand the pressure values Ap. and on the other hand 

the wheel radius values A^ obtained in step 504 and the wheel 
vibration values Af t obtained from the WVA 500 is as follows: 



AR 
A/ 




where 



(14a) 



-110 0 
-111-1 
0 0-11 



(15) 



is the modified version of the matrix H according to equation 
(2a) . The combined measurement vector y is determined from the 

values AJ^ and Af f in step 506, and the model matrix A is de- 
termined from the matrix H and the 4x4 identity matrix / 
which are both read from a matrix memory 510, from the parame- 
ter r Q which is the known nominal wheel radius read from the 
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nominal wheel radius memory 512, and from expansion coeffi- 
cients X r and X f read from a expansion coefficients memory 514. 

The model is based on a first order series expansion of the 
measurement values AR. and Af { in the tire pressure deviations 
Ap f / where the expansion coefficients X r and X f (which are 

preferably estimated in test drives) are the first order coef- 
ficients of this series expansion. Since typically a 3 0% tire 
pressure decrease approximately results in a 0.25% tire radius 
decrease, a reasonable choice of X r is: 

, -2.5 10" 3 

/L = — — , - (16) 

-0.3 

The choice of X f depends on the type of vibration phenomena 
measured by the WVA component . 

For those . skilled in the art, more complex relations between 
the measurement values , A/) and the tire pressure devia- 

tions Ap, can be contemplated. For this, a preferred general- 
ized functional dependence can be assumed as follows: 













g 2 (4p,,Ap 2 ,Ap 3 ,Ap 4 ) 


A ra 




£ 3 (Ap3,Ap 4 ) 


A/, 




/(Ap,) 


A/ 2 




/(Ap 2 ) 


A/ 3 




/(Ap 3 ) 


. A/ 4 _ 




/(Ap 4 ) 



where the arbitrary functions / , g x , g 2 and g 3 have to be es- 
timated by the contemplated model. 

In order to solve equation (14a) for the unknown tire pressure 
deviation values Ap, , the skilled person would contemplate 
various alternative solutions. 
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One preferred solution is based on a weighted least square es- 
timate which is performed in step 516 and derives weighted 
least square estimates Ap, for the pressure deviation vector 

5 Ap according to the following formula: 

JSp(i) = (A T A)' l A T y(t) . (18) 
Here, A is the matrix defined in equation (14a) . 

10 

In step 522, the pressure indication values r} { for each vehicle 
wheel ( i =1,2,3,4) can finally be calculated by assuming that 
the obtained pressure deviation values Ap. are Gaussian dis- 
tributed random quantities, which are computed in a preceding 
15 step 518 according to equation (10) with a standard deviation 
a Ap read from a standard deviation memory 520. Both steps 518 

and 522 commonly implement the following formula: 

77, =P/ r =l-0(-^-) . (19) 

20 

where <l>(x) is the Gaussian cumulative probability distribution 
function defined in equation (10) . The standard deviation value 
a Ap is preferably treated as a predefined tuning parameter 

which is obtained from test drives and most preferably stays 
25 constant during operation of the tire pressure estimating sys- 
tem. It could as well be determined from statistical computa- 
tions on the basis of subsequently derived pressure deviation 
values Ap f or by other ways which can be contemplated by those 
skilled in the art. 

30 

In step 524, tire pressure loss signals are then triggered in a 
similar way as in step 116 of the second embodiment by compari- 
son of the pressure indication values rj i with a critical 
threshold value. 



35 
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FIG. 6 shows a flow diagram of the process steps of an imple- 
mentation of the above described embodiment. In step 600 wheel 
radius values AR, are obtained from the WRA component . In a 
subsequent step 602, these wheel radius values Af^ are trans- 
formed to the modified wheel radius values A/^ / which are less 
sensitive to load changes. In step 604, vibration phenomena 
values A/, are obtained from the WA component. In step 606, 
the measurement vector y is constructed from the wheel radius 
values AR, and the vibration phenomena values hf. . The model 
matrix H , the nominal wheel radius r 0 and the coefficients X ff 
X r are fetched from a memory in step 608 to construct the model 
matrix A in step 610. Then, in step 612 the pressure parame- 
ters Ap. are estimated by a Weighted Least Square Method re- 
sulting in estimated pressure values Ap, . In step 614, the pre- 
determined standard deviation parameter a Lp is fetched from a 
standard deviation parameter memory and in step 616, the prob- 
ability value P* is computed from Ap, and a^ p according to 
equation (19) . In this embodiment, the pressure indication 
value ?] i is identical to the probability value P* (see equa- 
tion (19)) so that the calculation in step 618 is trivial. Fi- 
nally, in step 620, the pressure indication value 77., similarly 
to the first embodiment, is compared with a predetermined 
threshold value r\ c . If the pressure indication value 77, is less 
than the threshold value r\ c , a tire pressure loss signal is 
triggered in step 622. The procedure is repeated for every new 
data received from the WVA and WRA components. 

Besides, equation (14a) may be written in the following gener- 
alized, form which covers the general form of the WRA. measure 
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ment outputs according to equations (la-f) and other WRA meas- 
urement outputs as well as different WVA measurement outputs: 







~H 






I 



Ap 



(14b) 



where J denotes the identity matrix and H an arbitrary matrix 
adapted a specific WRA component (see above) . 



[Third embodiment] 



FIGS . 7 and 8 show a schematic diagram of a third embodiment of 
the invention and a flow diagram of the process steps of an im- 
plementation of the third embodiment, respectively. The third 
embodiment is based on the same model as the second embodiment. 
Instead of using a weighted least square estimate (step 516 of 
the second embodiment) for deriving the tire pressure deviation 
values Ap, from the vibration phenomena values Af f and the 

wheel radius values according to model equation (14a) (or 

wheel radius values AR { according to equation (14b) ) , alterna- 
tively an adaptive Kalman filter is used in a corresponding 
step 716. Therefore, a description of identical steps performed 
in the third embodiment is here omitted. Instead, reference is 
made to the description of the corresponding steps performed by 
the second embodiment as shown in FIGS. 5 and 6. 

The adaptive Kalman filter subsequently receives the input val- 
ues from the WVA and WRA components 700 and 702 (more precisely 
from the load change transformation step 704) and outputs its 
filter estimates for the pressure deviation values Ap, , whereby 
it. simultaneously adapts its filter parameters. The filter pa- 
rameter adaptation and the calculation of the filter estimates 
can be easily implemented within a pressure estimating system 
and offers less computational burden and memory demands . 

For a more detailed description of the adaptive Kalman filter 
see for example the books: 
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[1] Maybeck, Peter S.: Stochastic models, estimation, and con- 
trol, Mathematics in Science and Engineering, Volume 141, 1979; 
[2] Zarchan, Paul and Musoff, Howard: Fundamentals of Kalman 
Filtering: A Practical Approach, 2000, AIAA. 

The filter equation for the adaptive Kalman filter is given by 
equation (14a) according to: 

y(t) = AAp(t) + e(t) (2 0a) 

where e(t) is a noise vector. For the state sequence (i.e. the 
unknown tire pressure £p(t) ) a standard model is applied: 

Ap(* + 1) = Ap(t) + w(t) (20b) 

where w(t) is white noise. 

Applying models (20a) and (20b) , the Kalman filter then esti- 
mates a tire-pressure value Ap according to the following re- 
cursion : 

4ft0 = m - 1) -i- K(t)(y(t) - Abpit - 1)) 

K(t) = P(t - 1)A(s + A T P(t - 1)aY (20c) 
Pit) = P(t - 1) - P{t - l)A T (S + A T Pit - 1)^ 1 + Q 

The matrices S and Q are considered to be tuning parameters 
that control the adaptation speed of the Kalman filter. 

The pressure indication values rj i for each vehicle wheel 
(z =1,2,3,4) are then calculated from the obtained pressure de- 
viation values Ap, in the same way as shown in the second em- 
bodiment . 

FIG. 8 shows a flow diagram of the process steps of an imple- 
mentation of the above described embodiment. The process steps 
of the third embodiment differ from the process steps of the 
second embodiment described in FIG. 6 only in that the estima- 
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tion of the pressure parameters by a Weighted Least Square 
Method in step 612 of FIG. 6 is replaced by the above described 
Kalman Filter estimate in step 812 of FIG. 8. All remaining 
process steps are identical in both embodiments . 



[Fourth Embodiment] 



FIG. 9 shows a schematic diagram of a fourth embodiment of the 
invention. The fourth embodiment corresponds basically to the 
first embodiment with the exception that it adds prior to the 
wheel radius preprocessing step 106 of the first embodiment 
(referenced with the numeral 906 in FIG. 9) a load balance es- 
timation step 905. This load balance estimation step 905 is il- 
lustrated in more detail in FIG. 10. It is principally based on 
the model defined by equation (14b) . This model is modified in 
this embodiment by adding a load balance parameter / as a fur- 
ther component to the vector x which also comprises the cur- 
rent pressure values Ap f to be determined. The load balance pa- 
rameter / needs to be determined in order to compensate for 
load changes from the front wheels to the back wheels, and vice 
versa, during driving of the vehicle. The model matrix A read 
from the model matrix memory 1006 is modified with regard to 
the matrix A of equation (14a) in such a way, that the load 
balance parameter / appropriately influences the wheel radius 
deviation values and A^ . An increase of the value of / 

corresponds to a load shift to the back wheels whereas a de- 
crease of the value of / corresponds to a load shift to the 
front wheels. 
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It should be noted that in contrary to the second embodiment 
which is based on the modified tire radius values A/^ (see 
equation (14a) ) , the measurement vector y now includes the 
tire radius values A/?,, as obtained from the WRA component 1000 
output ting WRA measurement values according to equation (la) . 

In order to estimate the state parameters x (Ap. and /) in 
step 1014 , any one of the methods presented in the second and 
third embodiments can be applied (Least Mean Square estimation 
as performed in step 516, Kalman filter as applied in step 
716) . In particular, the Least Mean Square estimate of the load 
balance parameter / is given as follows: 

^^(A^r+A^ + ^A/.+A/.-A^-A/,) . (22) 

A load compensation of the value A/?, can then be obtained by 

subtracting the load balance estimation value / from the cor- 
responding tire radius values ARj in step 1016: 



AR = AR-IA 



(23) 



Then, further steps 906 to 916 which are identical to the steps 
106 to 116 of the first embodiment are performed. A description 
of theses identical steps is here omitted. Instead, reference 
is made to the description of the corresponding steps 106 to 
116 as shown in FIGS. 1 and 3. 

FIG. 11 shows a flow diagram of the process steps of an imple- 
mentation of the load balance estimation. In step 1100 and step 
1102, wheel radius values A/^. and vibration phenomena values 
A£ are obtained from the WRA component and the WVA component, 
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respectively. In step 1104, the matrix H , the nominal wheel 
radius r 0 and the coefficients X f , A, r are fetched from corre- 
sponding memories. In step 1106, the measurement vector y is 
constructed from the wheel radius values AR, and the vibration 
phenomena values Af. and, furthermore, the measurement vector y 
is normalized with the coefficients X f and X r . In step 1108, 
the model matrix A is constructed from the matrix H and the 
nominal wheel radius r 0 . In step 1110, the load and pressure 
deviation parameters / and Ap f are estimated, and by subtract- 
ing the estimated load value / from the wheel radius values 
AR, in step 1112 according to equation (23), the load corrected 

wheel radius values Ai^. are obtained, which are then output in 
step 1114. 

The remaining process steps of this embodiment are identical to 
the process steps of the first embodiment and it is for this 
referred to FIG. 3 and the corresponding description. 

The load change compensation does not need to be active at all 
time -instants. It may for example be triggered by a condition 

like A XL -& XR >1 & niin(|A J J,|A JK |)>J!r where k is a threshold. Using 
the above triggering condition, the functionality for single- 
tire deflations will not be affected at all. For axle-wise de- 
flations, the load compensation method will in principle re- 
quire that WVA as well as WRA indicate pressure loss before a 
"tire pressure loss" is issued. 

[Fifth Embodiment] 

In contrast to the fourth embodiment which corrects the impact 
of load changes with an estimated load parameter, the following 
alternative embodiment for a tire pressure estimation system 
treats these load changes as a random variable instead of an 
estimated variable. 
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Thus, similar to the "upper half" of equation (21) , AR may be 
related to Ap and the load balance parameter 1 according to 
the following equation by assuming that only Ap is a determi- 
nistic vector and 1 is a random variable: 



(A#//l r ) = — Hlsp + 



l + e 



(24) 



where e is a zero-mean Gaussian noise vector reflecting further 
uncorrelated statistical fluctuations of the tire radii. The 
10 noise introduced by the noise vector may be characterized by 

the following condition: E{ee T } = cr r 2 /» where a r 2 is the variance 

of the random variable e. In addition, it is assumed that the 
load balance 1 is a zero-mean Gaussian random variable which is 
not correlated with e. Thus, AR/A r is Gaussian distributed 
15 with mean H-hp/r^ and the following covariahce matrix Q m : 



(25) 



20 



25 



where erf is the variance of the random variable 2. 

The covariance matrix Qar ma Y alternatively be chosen as fol- 
lows, in particular for rendering the tire pressure estimation 
system more robust against drifts in the high precision yaw 
rate (cp. WO 01/87647) due to long curves: 



(26) 
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where af is the variance of the random variable 2. 

Accordingly, similar to the "lower half" of equation (21) , the 
normalized vibration phenomena vector Af is Gaussian dis- 

5 tributed with mean Ap and a covariance matrix Q^-cx^I • 

A (Maximum Likelihood) estimate for the pressure value Ap may 
then be obtained from the least squares solution of the follow- 
ing equation: 



10 



30 



es /a (Af/y). 



Ap (27) 



The estimated pressure value Ap may be further processed as 
described in the context of the second embodiment (equation 
15 (19)). Hereby, the parameters a?, a) and a] may be tuning pa- 
rameters or statistical parameters (see above) . 

[Sixth Embodiment] 

20 In the following, another alternative embodiment is described 
which relies on the standard x 2 ~ tes t . The basic assumption in 
this embodiment is that under the hypothesis "no deflation of 
tire i", Af i la f and Ar i la r are two zero-mean uncorrelated Gaus- 
sian random variables with unit variance, i.e. the sum of 

25 (Afi/o-j.) 2 and (Ar,/cr r ) 2 is % 2 -distributed with two degrees of 
freedom. 

Under this assumption, a simple x 2 - bas ed tire pressure estima- 
tion system may be based on the following relation: 




(28) 
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where, for each wheel i, A/) is the currently obtained value of 
the vibration phenomena , Ar, is the currently obtained value 
of the tire radius (see equation (6)), a f and a r are the cor- 
responding standard deviations, and k is a threshold value. 
5 Again, the standard deviations ay and a r may be predetermined 

tuning parameters obtained from test drives or "real" statisti- 
cal magnitudes which are actualized for each newly obtained 
value for the vibration phenomena and the tire radius. For 
those skilled in the art, various alternatives for choosing or 
10 measuring the standard deviation parameters are easily contem- 
plated, including those with dynamically adapted values. 

The threshold parameter k is chosen such that for example a 
flat tire condition will be detected if equation (28) is ful- 
15 filled. 

Besides, from £ for tire i (according to the left hand side of 
equation (28)) a pressure indication value rj { (similar to the 
one obtained in the above embodiments) can be obtained by the 
20 £ cumulated distribution function F(^,2) for two degrees of 
freedom 



77,=l-FCtr 2 ,2) (29) 

25 which is the confidence level for the hypothesis "no deflation 
of tire i". As an alternative of applying the above threshold 
value k together with equation (28) , the pressure indication 
value 77, obtained in equation (29) may then be treated simi- 
larly to the above embodiments. 

30 

[Seventh Embodiment] 

Below, an alternative embodiment for the tire pressure estima- 
tion system is described. The underlying method is non- 
35 parametric in the sense that no specific assumptions are intro- 
duced about the statistical distribution of Af and Ar . In- 
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stead, training data are collected during at least one test 
drive which, under normal drive conditions, are used to evalu- 
ate the "no deflation of tire i" hypothesis. The only assump- 
tion made is that the single tire pressure estimation for one 
single tire i can be performed by only considering Ar, and A/) 
corresponding to this tire i. 

FIG. 14 shows four different graphs, one for each of the four 
tires, in which ( A£ , Ar, ) parameter pairs obtained for each sin- 
gle tire i are plotted during a test drive under "no deflation 
of tire i" condition. A box plotted in the . graphs of FIG. 14 
indicates that this region in the ( Af { , Ar, ) parameter space has 
been "visited", which means that the corresponding ( Af. , A?;. ) pa- 
rameter pair occurred at least once during the recording of the 
training data. The areas in the graphs of FIG. 14 which are oc- 
cupied by boxes thus represent the "no deflation of tire i" 
condition. Consequently, a "deflation of tire i" condition 
shall only be indicated during normal drive condition for the 
case that a currently obtained ( bf. , Ar. ) parameter pair lies be- 
yond these areas filled with boxes. 

A simple tire pressure estimation system can be designed for 
example by specifying a preferably small number of grid points 
in the ( A/, r Ar f ) parameter space obtained for each tire i. These 
grid points are linked together and separate a first area rep- 
resenting the "deflation of tire i" condition from a second 
area indicating the "no deflation of tire i" condition. The 
tire pressure estimation system only has to decide if a cur- 
rently obtained ( Af i , Ar { ) parameter pair lies in the first or 
second area. This embodiment of a tire pressure estimation sys- 
tem is very cost-effective, in particular due to its low compu- 
tational complexity. A further advantage is that the process of 
separating the first and second areas (e.g. by using the above- 
mentioned grid points, etc.) may be easily automated. 
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[Eighth Embodiment] 

The following embodiment of a tire pressure estimation system 
is based on a different approach wherein the vibration phenom- 
ena and radius values from the WVA and WRA together with pres- 
sure deviation values are measured and collected during a num- 
ber of real -world driving situations (including cases with 
"deflation of tire i" and "no deflation of tire i" conditions) . 
For each time instant the following vector may be collected for 
the vibration phenomena and radius values from the WVA and WRA: 

Furthermore, a general model relating these vectors <p(t) with 

the pressure deviation values Ap(f) may be established as fol- 
lows : 

4p(0 = /<*(0.*) < 31 > 

where f is an arbitrary function and the generalized variable 0 
designates some unknown parameters. One preferable choice for 
the function f is a series expansion in the vector <p(t) with 

expansion parameters 9. For example neural networks, radial ba- 
sis function networks, and so forth may be used as such a se- 
ries expansion. 

The unknown parameters 6 may then be estimated on the basis of 
a given function and a collection of N measurements by minimiz- 
ing the following standard least -squares criterion: 

T 

0 = argmin, V N (0) = aigmm,~2[4Pto^ (32) 

where the vectors <p(t) and the pressure deviation values Ap(t) 
are measured data values, and N designates the number of such 
measured data values. A decision between a "deflation of tire 
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i" condition and a "no deflation of tire i" condition is fi- 
nally based on the outputs of the estimated model : 

4K0 = /(rt0A (33) 

For the measured pressure deviation values Ap(f) / one may use 
binary values to represent pressure deviations from the nominal 
pressure, e.g., let Ap.(*) = 0 represent a non significant pres- 
sure drop situation and Ap { (t) = l a significant pressure drop 

situation in tire i. These binary values may be manually en- 
tered into the system by driving under a "no deflation of tire 
i" condition (Ap t (t) = 0 will be entered for all measured vectors 
<p(t)) and under a "deflation of tire i" condition (Ap.(0 = l will 
be entered for all measured vectors <p(f)) . Alternatively, fixed 
pressure deviation values Ap(f) may be entered where the pres- 
sure is changed in a step-like manner. As another alternative, 
a tire pressure measuring system may be used to enter the meas- 
ured pressure deviation values A?(0 i nto t!h& system. 

[Miscellaneous] 

According to an exemplified embodiment of the invention, the 
tire pressure estimating system may be any machine capable of 
executing a sequence of instructions that specify actions to be 
taken by that machine for causing the machine to perform any 
one of the methodologies discussed above. The machine may be an 
application specified integrated circuit (ASIC) including a 
processor and a memory unit- The instructions may reside, com- 
pletely or at least partially, within the memory and/or within 
the processor. 

In particular, the tire pressure estimating system may be im- 
plemented in the form of a computer system within which a se- 
quence of instructions may be executed. The computer system may 
then further include a video display unit, an alpha-numeric in- 
put device (e.g. a keyboard), a cursor control device (e.g. a 
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mouse) , a disk drive unit. The disk drive unit includes a ma- 
chine-readable medium on which is stored the sequence of in- 
structions (i.e., a computer program or software) embodying any 
one, or all, of the methodologies described above. 

5 

The computer program product may be a machine-readable medium 
which is capable of storing or encoding the sequence of in- 
structions for execution by the machine and that cause the ma- 
chine to perform any one of the methodologies of the present 
10 invention. The machine-readable medium shall accordingly be 
taken to include, but not be limited to, solid-state memories, 
optical and magnetic disks, and carrier wave signals. 

In the above description, for simplification, it is mostly re- 
15 f erred to the method case. The system case can be easily de- 
rived from the latter one by replacing at the appropriate 
places in the description the expression f step f for the method 
case by the expression f component 1 for the system case. 

20 Furthermore, all other publications and existing systems men- 
tioned in this specification are herein incorporated by refer- 
ence. 

Above, various embodiments of the inventions are described 
25 within the context of vehicle wheels. However, it is noted that 
the invention is not limited to vehicle wheels but may include 
other wheels, as airplane wheels, transport wheels, machine 
wheels and the like. Furthermore, it is not limited to a par- 
ticular number of wheels. 

30 

The tire pressure estimation system may further comprise detec- 
tion means for detecting driving on gravel. On gravel the sys- 
tem may be tuned such that it is less sensitive to pressure 
changes . In principle, only "alarms on a single tire may be al- 
35 lowed on gravel. 
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It is noted that the WRA and WA components are not limited to 
the specific WRA and WVA components as described above. 

In general, any WRA component may be used which provides single 
5 wheel radius values or (linear) combinations thereof. For exam- 
ple , single wheel radius values or combinations thereof may be 
provided by combining the absolute vehicle velocity with the 
rotational velocity of a tire. The absolute vehicle velocity 
may for example be obtained 
10 - by correlating wheel speed signals of front and rear 
wheels (for more details reference is made to the co- 
pending application "Determination of a vehicle's absolute 
velocity" of the same applicant, the content of which is 
herein incorporated) , 
15 - by correlating road surface pictures captured by cameras 
near the front and rear tires, 
from navigation system data, 
from a "fifth" (e.g. non-def ormable) * wheel, 
from observations of road markings, 
20 - from detection of a magnetic tape on the road surface, or 
from lateral and longitudinal dynamics, etc. 

More precisely, in the first mentioned velocity determination 
of a vehicle having at least one pair of a front and a rear 

25 wheel which are spaced by a wheel spacing jB, front and rear 
wheel speed signals are determined which are indicative of the 
time dependent behavior of the front and rear wheel speeds, re- 
spectively. The front and rear wheel speed signals are corre- 
lated in order to determine a specific correlation feature in- 

30 dicative of the time delay between the front wheel and rear 
wheel speed signals. The velocity of the vehicle is determined 
based on said correlation feature and the wheel spacing. In 
particular, the action of correlating is based on a correlation 
function of the front and rear wheel speed signals, wherein 

35 said correlation function is a function of a time difference, 
and the correlation feature is the specific time difference 
which corresponds to the maximum of said correlation function. 
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Thereby, the vehicle velocity is computed from the specific 
time difference. For example, the correlation function i^ 3 (r) 
for the front-left and the rear left wheel may be defined as 

(*) = *>i (0 - Efa (0)X«3 i*-r)- Efa 0 ~ *)))] • ( 3 4 > 

with denoting the wheel speed signal of the respective 

wheel i, £(©(0) denoting the expectation value of a>(t) , and r 
denoting the time delay. Small disturbances injected by uneven 
road surface will occur first on o>,(f) and B/v t seconds later on 
a> 3 (t) . Here v, denotes the velocity of the car. The correlation 
function i^ 3 (r) is a function of time delay % and will show a 
peak at 

r = — (35) 
v, 

When driving with constant speed, one can directly obtain 
J^ s (r) . The time delay is obtained from -R, 3 (V) by t = argmax i?, 3 (r) . 

Inserting the obtained i in equation (35) then yields the ab- 
solute velocity of the car. 

Alternatively, the action of correlating comprises a compensa- 
tion with a wheel speed signal of the front or the rear wheel 
such that the correlation feature remains unchanged with vary- 
ing vehicle velocity but changes with varying wheel radius of 
the respective wheel. Hereby, a correlation function of the 
front and rear wheel speed signals is applied which is a func- 
tion of the reciprocal of the product of the known wheel speed 
signal and the respective unknown wheel radius. Then, the cor- 
relation feature is the specific wheel radius which corresponds 
to the maximum of the correlation function, whereby the vehicle 
velocity is computed from the specific wheel radius and the 
corresponding wheel speed signal. In detail, using the relation 
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r-B/co^ and defining x = 2?/r, in equation (34) results in a ve- 
locity compensated correlation function 

j? 13 (x) = Efa (0 - Efa (0)X«s (' - *M (0) - JE(®3 0 - *M (0))) • < 3 6 ) 

By maximizing J^ 3 (x) with respect to the variable x = B/r l one can 
obtain an estimation value r x for the wheel radius of the 

JD 

front-left wheel £ = argmax i? 13 (— ) . Analogously, this allows to 

* r 

estimate all wheel radii r t of the vehicle even at rapidly 
varying vehicle velocities. From these wheel radii r t , the ab- 
solute velocity v. =<B.r. at each wheel can be determined and sim- 
ple geometrical transformations can be used to find the veloc- 
ity at any position of the vehicle. Naturally, the wheel radii 
r t thus determined may be directly used for a combination with 
the wheel vibration data values. 

Alternatively, the correlation feature is obtained by 
Fourier transforming the wheel speed signals to obtain Fourier 
transformed wheel speed signals and calculating a phase func- 
tion of the ratio of the Fourier transformed wheel speed sig- 
nals of the front and the rear wheel. Thereby, the slope of 
said phase function is the correlation feature indicative of 
the time delay. In detail, Fourier transforming the wheel speed 
signals co^t) with standard FFT-methods results in corresponding 
Fourier coefficients A time shift of r of the wheel 

speed signals in the time domain a>(f-r) corresponds to a 

phase shift £!(/) -» Q{f)e~ anft in the frequency domain. One ap- 
proach is to adjust a straight line to the phase curve of 
which should have a slope 2ttt . 

Similarly, the WVA component could be any component which pro- 
vides wheel vibration data values, for example by using 
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microphones for recording driving noise near the vehicle 
wheels, 

accelerometers for measuring vehicle vibrations, or 
wheel suspension sensors, etc. 

Other input signals can also be used in the tire pressure esti- 
mating system for load compensation. Any signal that measures 
the weight on a wheel or the entire vehicle may be contem- 
plated, such as wheel suspension, trailer-sensor, estimator of 
physical mass etc. 

Although certain methods, systems and products constructed in 
accordance with the teachings of the invention have been de- 
scribed herein, the scope of coverage of this patent is not 
limited thereto. On the contrary, this patent covers all em- 
bodiments of the teachings of the invention fairly falling 
within the scope of the appended claims either literally or un- 
der the doctrine of equivalents. 



